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ABSTRACT 

A new analytical solution is introduced for the effect of chemical reaction on mixed convective 
heat and mass transfer in a vertical double passage channel. The vertical channel is divided into two 
passages (by means of a baffle) for two separate flow streams. Each stream has its own individual 
velocity, temperature and concentration fields. After placing the baffle the fluid is concentrated in one 
of the passage. Approximate analytical solutions are found for the coupled nonlinear ordinary 
differential equations using regular perturbation method (PM) and Differential Transform method 
(DTM). The validity of the Differential Transform series solutions are verified with the regular 
perturbation method. The velocity, temperature and concentration solutions are obtained and 
discussed for various physical parameters such as thermal Grashoff number, mass Grashoff number, 
Brinkman number and chemical reaction parameter at different positions of the baffle. It is found that 
the thermal Grashoff number, mass Grashoff number, Brinkman number enhances the flow whereas 
chemical reaction parameter reduces the flow at all baffle positions. It is also found that as Brinkman 
number increases the DTM and PM show more error. 

KEYWORDS: Baffle, first order chemical reaction, mixed convection, perturbation method, 
Differential Transform method. 



I. INTRODUCTION 

Study of mixed convection in the channel has been to the focus of lot of investigation during the last 
three decades because of the multiple applications in which it is involved. These includes cooling of electronic 
equipment, heat exchangers, chemical processing equipment, gas-cooled nuclear reactors and others. Tao [1] 
analyzed the laminar fully developed mixed convection flow in a vertical parallel -plate channel with uniform 
wall temperatures. Aung and Worku [2, 3] discussed the theory of combined free and forced convection in a 
vertical channel with flow reversal conditions for both developing and fully developed flows. The case of 
developing mixed convection flow in ducts with asymmetric wall heat fluxes was analyzed by the same authors 
[4]. Recently, Prathap Kumar et al. [5] and Umavathi et al. [6, 7] studied the mixed convective flow and heat 
transfer in a vertical channel for immiscible viscous fluids. 

The rate of heat transfer in a vertical channel could be enhanced by using special inserts. Heat transfer 
in such partially divided enclosures has received attention previously due to its applications to design energy 
efficient buildings and reduction of heat loss from flat plate solar collectors. When the channel is divided into 
several passages by means of plane baffles, as usually occurs in heat exchangers or electronic equipment, it is 
quite possible to enhance the heat transfer performance between the walls and fluid by the adjustments of each 
baffle position and strengths of the separate flow streams. In such configurations, perfectly conductive and thin 
baffles may be used to avoid significant increase of the transverse thermal resistance. For a number of fluids, the 
density-temperature relation exhibits an extreme. Because the coefficient of thermal expansion changes signs at 
this extremum. Simple linear relations for density as a function of temperature are inadequate near the 
extremum. Dutta and Dutta [8] first reported the enhancement of heat transfer with inclined solid and perforated 
baffles. Later Dutta and Hossian [9] did the experimental study to analyze the local heat transfer characteristics 
in a rectangular channel with inclined solid and perforated baffles. Salah El-Din [10, 11] published a series of 
papers on mixed convection in a vertical channel by introducing a perfectly conducting baffle. 
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Mousavi and Hooman [12] studied numerically the fluid flow and heat transfer in the entrance region 
of a two dimensional horizontal channel with isothermal walls and with staggered baffles. Heat transfer 
enhancement in a heat exchanger tube by installing a baffle was reported by Nasiruddin and Siddiqui [13]. They 
found that the average Nusselt number for the two baffles case is 20% higher than the one baffle case and 82% 
higher than the no baffle case. Recently, Prathap Kumar et al. [14, 15] studied the flow characteristics of fully 
developed free convection flow of a Walters fluid (Model B') in a vertical channel divided into two passages. 
Umavathi [16] analyzed the effect of the presence of a thin perfectly conductive baffle on the fully developed 
laminar mixed convection in a vertical channel containing micropolar fluid. 

Combining heat and mass transfer problems with a chemical reaction are of importance in many 
processes and have, therefore, received a considerable amount of attention in recent years. In such processes as 
drying, energy transfer in a wet cooling tower, and the flow in a desert cooler, heat and mass transfer occurs 
simultaneously. Mixed convection processes involving the combined mechanisms are also encountered in many 
natural processes, such as evaporation, condensation, and agricultural drying, and in many industrial 
applications, such as the curing of plastics and the manufacture of pulp-insulated cables [17]. In many chemical 
engineering processes, chemical reactions take place between a foreign mass and the working fluid which 
moves due to the stretch of a surface. 

The order of the chemical reactions depends on several factors. One of the simplest chemical reactions 
is the first-order reaction in which the rate of the reaction is directly proportional to the species concentration. 
Chamkha [18] studied the analytical solutions for heat and mass transfer by the laminar flow of a Newtonian, 
viscous, electrically conducting and heat generating/absorbing fluid on a continuously moving vertical 
permeable surface in the presence of a magnetic field and the first-order chemical reaction. Muthucumaraswamy 
and Ganesan [19] studied the numerical solution for the transient natural convection flow of an incompressible 
viscous fluid past an impulsively started semi-infinite isothermal vertical plate with the mass diffusion, taking 
into account a homogeneous chemical reaction of the first order. 

The coupled nonlinear ordinary differential equations governing the flow are solved using regular 
perturbation method which is the oldest method used by many researchers. In this paper a new method known as 
Differential Transform method is applied to find the analytical solution. The main advantage of DTM is that it 
can be applied directly to nonlinear differential equations without requiring linearization, discritization, or 
perturbation. This method is well addressed in [20-24]. Recently Umavathi et al. [25] solved the coupled 
nonlinear equations governing the flow for magnetoconvection in a vertical channel for open and short circuits 
usng Differential Transform method. The aim of this paper is to investigate effect of first order chemical 
reaction of viscous fluid in a vertical channel in the presence of a thin conducting baffle. After inserting the 
baffle, the fluid in stream- 1 is concentrated. Analytical solutions are found using PM and using DTM. 

II. MATHEMATICAL FORMULATION 

Consider a steady, two-dimensional laminar fully developed free convection flow in an open ended 
vertical channel filled with purely viscous fluid. The X-axis is taken vertically upward, and parallel to the 
direction of buoyancy, and the 7-axis is normal to it. Walls are maintained at a constant temperature and the 
fluid properties are assumed to be constant. The channel is divided into two passages by means of thin, perfectly 
conducting plane baffle and each stream will have its own pressure gradient and hence the velocity will be 
individual in each stream. 
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Figure 1. Physical configuration. 



The governing equations for velocity, temperature and concentrations are 
Stream-I 

dP d 2 U. 





, dP d'U. 

PSP T {T i -T w )+ pgp c {C x -C^ +M i 

J/ dX dY 



d_\ _v_ f dU ^ 





dY z aC A dY 
d 2 C 

D kC = 

dY 2 
Stream-II 

x dP d 2 U , 

P80 T { T i- T w t )- +^ f 

v dX dY 

d 2 T 2 v ( dU 2 \ 

T+ 1 ~\ = 

dY aC p { dY ) 



(1) 

(2) 
(3) 

(4) 
(5) 



subject to the boundary and interface conditions on velocity, temperature and concentration as 

U = ,T, = T„, , C = C, , at Y = -h 



U, = 0,r 2 = T w , at Y 



U = ,U = 0, T = T ,- 



,C = C 2 ,at Y = h 



(6) 



Introducing the following non-dimensional variables, 



U 



U , 



T - T„ 



T -T 

w, 



gPr^Th gP c ^Ch : 

'' G r = 2 ' G C = 2 



C, - C„ u 



P = — = 

//[/, dX 

where ( = 1,2. 
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AT = T - T AC = C - C Y* = Y = — 



Br 
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kAT 



(7) 



Issn 2250-3005 



|October| I2013 



Page 28 



Perturbation Technique And Differential. 



The momentum, energy and concentration equations corresponding to stream-I and stream-II become 
Stream-I 



d 2 u. 



dy- 



GR r 0, + GR 







d 2 0, ( du, } 
+ Br 



dy' 



dy' 
Stream-II 



dy 



a V = 



d 2 w„ 



dy 



, + GR T 2 - p = 



d 6> 



fir 



f fif« 2 1 
dy 



(8) 

(9) 
(10) 

(11) 
(12) 



subject to the boundary conditions, 

u x = , t = 1 , </> = 1 ,at y = -1 

, at y = 1 



«j = , 



w j — , k 2 = , 6 1 , = 6» 2 , 



where a 



D 



c, - c„ 



n ,at y 



(13) 



III. SOLUTIONS 

The exact solution for concentration distribution is found using Eq. (10) and is given by 

B^Cosh (a y) + B 2 Sinh (a y) 



(14) 



3.1 Perturbation Method 

Equations (8), (9), (11) and (12) are coupled non -linear ordinary differential equations. Approximate 
solutions can be found by using the regular perturbation method and Differential Transform method. The 
perturbation parameter is considered as Brinkman number B r . Adopting this method, solutions for velocity and 
temperature are assumed in the form 

u i(y) = u w (y) + Bru n (y) + Sr2 ", 2 (y) + ■•• (1 5 ) 



9 ,{y) = 6 >w(y) + Br0 n{y) + Brl0 i2 ( y) + ■■■ 

where the subscript i = 1 and 2 represents stream-I and stream-II respectively. 



(16) 



Substituting Eqs. (15) and (16) into Eqs. (8), (9), (11) and (12) and equating the coefficients of like 
power of B r to zero and one, we obtain the zeroth and first order equations as 
Stream-I 

Zeroth order equations 



d 2 0, 



dy' 

d 2 u, 







dy 



2 , GR T lB 



GR 



(17) 
(18) 
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First order equations 

d 2 9, , ( du,.^i 



dy 2 I dy 



d Mj 

77 



dy 1 



First order equations 



(19) 



-+ GR T n = (20) 



Stream-II 

Zeroth order equations 

— =0 (21) 

dy 

d 2 u la 

+GR T 2O - p = (22) 



d 2 2l ( du 2n ~\ 

2J -+ I — 2JL I = (23) 



dy' { dy 
d'u 

^+GR T 21 = O (24) 

dy 

The corresponding zeroth order boundary conditions reduces to 
«io = °. ^„ =l,at y = -1 
«20 = , # 20 = , at y = 1 

«io = °. "20 =0, f lfl = # 20 , -^=^,aty=>.* (25) 

iiy cfy 

The corresponding first order boundary conditions reduces to 
u u = , <9 U = at y = -1 

Mji = , 21 = at y = 1 

b u = 0, u 2l = , U = 21 , -= ^aty=y* (26) 

The solutions of zeroth and first order equations (17) to (24) using the boundary conditions as given in 
Eqs. (25) and (26) are 

Zeroth-order solutions 
Stream-I 

°io = c ,y + c i (27) 

u io = A 2 + A { y + r^y~ + r 2 y y + r 4 C osh (a y) + r 5 Sinh (a y) (28) 

Stream-II 

6 20 = + C 4 (29) 

"20 = A 4 + \y + r ^y 2 + r y (30) 
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First order solutions 
Stream-I 

e n = G z + G iy + p^y 2 + p 2 y 3 + p^y* + p^ + p 5 y b + p b c° sh i 2a >') 

+ p 7 Sinh ( 2 a y ) + p g Cosh ( a y ) + p 9 Sinh (a y ) + p 10 yC osh ( a _y ) (31) 
+ p u yS;n/i(ay) + Ppj'Cos^a j) + p B y ! 5 in h (a y ) 
« u = G 6 + G s y + + R 2 y' + R,y' + R 4 y 5 + fl 5 y 6 + R 6 y 7 + R y y* 

+ R g Cosh (2a y) + R g Sinh (2a y) + R l0 Cosh (a y ) + R n Sinh (a y ) (32) 
+ R lz yCosh (a y ) + /? 13 yS inh (a y ) + R l4 y 'Cosh (a y ) + /? 15 y z Sinh (a y ) 

Stream-II 

2i = G 4 + G 3? + ?,:y 2 + ? 2 y 3 + q 3 y 4 + q 4 y 5 + q 5 y" (33) 

« 21 = G 8 + G 7 y + 5,y 2 + S 2 / + 5 3 y 4 + 5 4 y 5 + S 5 y " + S b y 1 + S 7 y* (34) 

3.2 Basic concepts of the differential transform method 

The analytical solutions obtained in Section 3.1 are valid only for small values of Brinkman number 
B r . In many practical problems mentioned earlier, the values of Sc are usually large. In that case analytical 
solutions are difficult, and hence we resort to semi-numerical-analytical method known as Differential 
Transform method (DTM). The general concept of DTM is explained here: The k lh differential transformation of 

an analytical function F ( k ) is defined as (Zhou [20]) 

*"(*) = —! )r-\ ' < 35 ) 

k\\ drj 

and the inverse differential transformation is given by 

f(l) = Y J F(k)(r 1 - rlo )\ (36) 

k =0 

Combining Eqs. (35) and (36), we obtain 
kl di] k 

From Eqs. (35)— (37), it can be seen that the differential transformation method is derived from Taylor's 
series expansion. In real applications the sum V F ( & ) (?7 - t] ) is very small and can be neglected when n 

k=n 

is sufficiently large. So / (77 ) can be expressed by a finite series, and Eqn. (36) may be written as 



/(*) = ! 



(37) 



/(*) = E F ( k )(n-ri ) , (38) 

k=0 

where the value of n depends on the convergence requirement in real applications and F (k ) is the 

differential transform of / (77 ) . Table 1 lists the basic mathematics operations frequently used in the following 
analysis. 
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Table 1 The operations for the one-dimensional differential transform method. 



Original function 


Transformed function 


yin) = g(n) ± h (n) 


Y (k) = G(k) ± H (k) 


y(rj) = a g(rj) 


Y(k) = a G(k) 


dg (r] ) 

y(n) - 

d r] 


Y(k) = (k + l)G(k + 1) 


y(n)- 2 

d r] 


Y(k) = (k + l)(k + 2)G(k + 2) 


yin) = g(i)h(n) 


k 

Y(k) = ^ G(l)H (k - I) 

1=0 


yi.v) = v"' 


[ 1, if k = m 
Y (k) = 5 (k - m) = { 

[0, if k m 



Taking differential transform of Eqs. (8), (9), (11) and (12), one can obtain the transformed equations as 
Stream-I 



©,(* + 2) 
d> (k + 2) = 



1 



(* + l)(k + 2) 
Br 

(k + l)(k + 2) \ 
a 2 <t> (k) 

+ + 2) 



(GR T @ l (k) + GR c O (k) - p 8 (k)) 

k 

£ (k - r + l)(r + l)U 1 (k - r + 1 ) C/ ± (r + 1) 



(39) 
(40) 
(41) 



Stream-II 

t7 2 (Jfc + 2) 

© 2 (Jfc + 2) 



1 



Jfc + + 2) 

Br 



GR T ® z (k) - pS(k)) 

k 

^ (k - r + l)(r + 1)U , (k - r + 1)U 2 (r + 1) 



(42) 
(43) 



(Jfc + l)(k + 2) ^ 

where, U 1 (k) , U 2 ( fe ) , 2 ( £ ) and <1> ( k ) are the transformed notations of u x ( y ) , « ( y ) , 

r i, if jfc = o 

( - v ) ' ^2 ( >' ) an d ( y ) respectively. 5 (k ) = <{ 

[0, if k > 

The following are the transformed initial conditions 

^ (°) = c , >u t (1) = c 2 , f7 2 (0) = c 3 ,U 2 (1) = C4 , 

©,(0) = d lt 0,(1) = d 2) 8,(0) = rf 3 , 2 (1) = rf 4 , 

O (0) = e, , © (1) = e 2 (44) 



Using the boundary condition (13), we can evaluate c, , c 2 , c , c 4 , d x , d , , rf , and e 2 . 
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IV. RESULTS AND DISCUSSIONS 

The objective of the present study is to understand the characteristics of mixed convection of a viscous 
fluid in a vertical double passage channel in the presence of chemical reaction. The solutions are found using 
perturbation method and Differential Transformation method. The physical parameters such thermal Grashoff 

number GR T , mass Grashoff number GR C , Brinkman number Br (or perturbation parameter) and chemical 

reaction parameter a , are fixed as 5, 5, 0.1, and 0.5 respectively, for all the graphs except the varying one. The 
effect of these parameters on velocity, temperature and concentration are shown in Figs. 2-10. The effect of 

thermal Grashoff number GR T (ratio of Grashoff number to Reynolds number) on the velocity and temperature 

is shown in Figs. 2a,b,c and Figs. 3a,b,c at all three different baffle positions (i.e. y* = -0.8 , 0.0 and 0.8). As 
the thermal Grashoff number increases, the velocity and temperature increases at all the baffle position whereas 
the maximum velocity field is observed in the wider stream. It is also observed form Figs. 3a,b,c that the 
temperature distribution is more effective near the left wall when compared to right wall. Further it is well- 
known that since Grashoff number is the ratio of buoyancy force to viscous force, increase in Grashoff number 
is to increase the buoyancy force and hence increases the concentration also. Therefore as the thermal 
Grashoff number increases velocity and temperature increases at all baffle position in both the streams. The 

effect of mass Gerashof number GR C (ratio of modified Grashoff number to Reynolds number) is shown in 
Figs. 4a,b,c for velocity field and in Figs. 5a,b,c for the temperature field. Here also the effect of GR C is to 
increase the velocity and temperature field in both the streams. It is seen from Figs. 4a and 5a (y* = -0.8 ) that 
the effect of GR C on the velocity and temperature fields is not effective whereas when the baffle position is at 

y* = 0.0 and 0.8 the flow field is enhanced as GR C increases. The similar result is also observed by 
Fasogbon [26] for irregular channel. 

The effect of Brinkman number Br on the velocity and temperature fields are shown in Figs. 6a,b,c 
and Figs. 7a,b,c respectively. As the Brinkman number increases, both the velocity and temperature increases in 
both the streams at all baffle positions. One can see from temperature equation that increase in Brinkman 
number increases the viscous dissipation and hence the temperature increases, which intern influences the 
velocity and temperature. The effect of first order chemical reaction parameter a , on the velocity, 
temperature and concentration fields is shown in Figs. 8a,b,c, Figs. 9a,b,c and Figs. 10a,b,c respectively. As a 
increases the velocity and temperature decreases in stream-I, and remains invariant in stream-II when the baffle 
position y * = -0.8 . But when the baffle position is at y * = & . 8 the effect of a is more effective in 
stream -I and less effective in stream -II. This is because the fluid is concentrated in stream-I only. The effect of 
chemical reaction parameter a is to decrease the concentration distribution as seen in Figs. 10a,b,c, which is 
the similar result obtained by Srinivas and Muturajan [27] for mixed convective flow in a vertical channel. It is 
observed from Tables 2a, 3a and 4a that results of DTM and PM agree well in the absence of Brinkman number 

at all the baffle positions. For large values of Brinkman number [Br ^ ) , DTM and PM solutions show 

difference as seen in Tables 2(b,c) to 4(b,c). It is also observed from these tables that the error of DTM and PM 
is very less in smaller stream when compared to bigger stream at all baffle position for B r * . 

V. CONCLUSION 

The effect of first order chemical reaction in a vertical double passage channel filled with purely 
viscous fluid was investigated. The solutions of the governing equations and the associated boundary conditions 
have been obtained by using regular perturbation method and differential transform method. Main findings are 
summarized as follows: 

[1] Increasing thermal Grashoff number, mass Grashoff number and Brinkman number increases the velocity 

and temperature in both the streams at all different baffle position. 
[2] Increase in the chemical reaction parameter suppresses the velocity and temperate in stream-I and remains 

invariant in stream-II. 

[3] The use of baffle in the flow channel resulted in the heat transfer enhancement as high as compared to the 

heat transfer in a channel without baffle. 
[4] Chemical reaction parameter was to decrease the flow field. 

[5] An excellent agreement was observed with the results of DTM and PM for small values of Brinkman 
number. 
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NOMENCLATURE 



Br 

Ci 
C o 

c 

p 

c 

p 
D 



Brinkman number 



( 2 
1 "i M, 1 



K AT 

v J 



the Concentration in Stream-I 

reference concentration 

specific heat at constant pressure 

dimensionless specific heat at constant pressure 

diffusion coefficients 
acceleration due to gravity 
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( h'g PAT } 
G r Grashoff number I I 

I v 1 J 

( g fic AC h' ^ 

G c modified Grashoff Number | I 

I r 1 ) 

GR T thermal Grashoff number (= Gr I Re ) 
GR C mass Grashof number (= Gc I Re ) 
h channel width 

h width of passage 

k thermal conductivity of fluid 

h dp 

p non-dimensional Pressure Gradient (= ) 

U lM dX 

U h 

Re Reynolds number ( ) 

r 

T , T 2 dimensional temperature distributions 
r , r temperatures of the boundaries 

U , reference velocity 

U 1 ,U dimensional velocity distributions 

u. , u. non dimensional Velocities in Stream-I, Stream-II 

y baffle position 

GREEK SYMBOLS 

a chemical reaction parameters 

P coefficients of thermal expansion 
P c coefficients of concentration expansion 
AT , AC difference in Temperatures & Concentration 

£ perturbation Parameter 

( T.-T. } 

6 non-dimensional temperature 

Y kinematics viscosity 

</> non-dimensional concentrations 

p density 

p viscosity 

SUBSCRIPTS 

i refer quantities for the fluids in stream-I and stream-H, respectively. 
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Fig. 2: velocity distribution for different values of thermal Grashof number C i R ^ 

at Ca) y*=-0.8 (b)v*=0.0 Cc) v *=0.8 
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Figu.3: Temperature profile for different values of ratio of Grashof number to 
Reynolds number GR at (a)y*=0.8 (b)y*=0 (c)y*=0.8 
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Kig.4: Velocity profile for different values ofratio of modified Grashoff number to 
Reynolds number GR„ at (a)y*=-0.8 (b)y*=0 (c)y*=0.8 
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5: Temperature profile for different values of ratio of modified Grashof number to 
Reynolds number GR at (a)y*=-0.8 (b)y*=0 (c)y*=0.8 
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Fig. 8: Velocity profile for different values of chemical reaction parameter ot 

at (a) v *=-0.8 CbOv*=0 (c) v*=0.8 
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Fig. 9: Temperature profile for different values of chemical reaction parameter ot 
at (a) y*=-0.8 (b)y*=0 (c) y*=0.8 
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Figure lO. Concentration profile for different values of chemical reaction parameter ot 



Table 2a Comparison of velocity and temperature with Br = o , GR T = 5 , GR C = 5 , p = -5 and 



0.0 . 





Velocity 


Temperature 


y 


DTM 


PM 


EiTor 


DTM 


PM 


Error 


-i 








0.0000 


1.000000 


1.000000 


0.0000 


-0.75 


1.266461 


1.266461 


0.0000 


0.875000 


0.875000 


0.0000 


-0.5 


1.659656 


1.659656 


0.0000 


0.750000 


0.750000 


0.0000 


-0.25 


1.227398 


1.227398 


0.0000 


0.625000 


0.625000 


0.0000 











0.0000 


0.500000 


0.500000 


0.0000 


0.25 


0.605469 


0.605469 


0.0000 


0.375000 


0.375000 


0.0000 


0.5 


0.781250 


0.781250 


0.0000 


0.250000 


0.250000 


0.0000 


0.75 


0.566406 


0.566406 


0.0000 


0.125000 


0.125000 


0.0000 


1 








0.0000 








0.0000 



Table 2b Comparison of velocity and temperature with 



Br =0.05 , Gff T 



-5 and 



0.0 





Velocity 


Temperature 


y 


DTM 


PM 


Error 


DTM 


PM 


Error 


-i 








0.0000 


1.000000 


1.000000 


0.0000 


-0.75 


1.339968 


1.329565 


0.0104 


0.989529 


0.973754 


0.0158 


-0.5 


1.771965 


1.755951 


0.0160 


0.933166 


0.907116 


0.0261 


-0.25 


1.321337 


1.307845 


0.0135 


0.870308 


0.834778 


0.0355 











0.0000 


0.761836 


0.722594 


0.0392 


0.25 


0.682521 


0.670711 


0.0118 


0.583393 


0.551573 


0.0318 


0.5 


0.870491 


0.856765 


0.0137 


0.393647 


0.371510 


0.0221 


0.75 


0.622964 


0.614236 


0.0087 


0.202259 


0.190149 


0.0121 


1 








0.0000 








0.0000 
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Table 2c Comparison of velocity and temperature with Br = 0.15 , GR T = 5 , GR C = 5 , p = -5 and y* = o.o . 





Velocity 


Temperature 


y 


DTM 


PM 


Error 


DTM 


PM 


Error 


-i 








0.0000 


1.000000 


1.000000 


0.0000 


-0.75 


1.651154 


1.455775 


0.1954 


1.465379 


1.171263 


0.2941 


-0.5 


2.249870 


1.948541 


0.3013 


1.711319 


1.221349 


0.4900 


-0.25 


1.723060 


1.468738 


0.2543 


1.925328 


1.254333 


0.6710 











0.0000 


1.915429 


1.167783 


0.7476 


0.25 


1.027983 


0.801196 


0.2268 


1.514892 


0.904720 


0.6102 


0.5 


1.271710 


1 .007795 


0.2639 


1.040516 


0.614531 


0.4260 


0.75 


0.877917 


0.709895 


0.1680 


0.554797 


0.320447 


0.2344 


1 








0.0000 








0.0000 



Table 3a Comparison of velocity and temperature with Br = o , GR T = 5 , GR C = 5 , p = -5 and y* = -0.8 . 





Velocity 


Temperature 


y 


DTM 


PM 


EiTor 


DTM 


PM 


Error 


-i 








0.0000 


1.000000 


1.000000 


0.0000 


-0.95 


0.055395 


0.055395 


0.0000 


0.975000 


0.975000 


0.0000 


-0.9 


0.073646 


0.073646 


0.0000 


0.950000 


0.950000 


0.0000 


-0.85 


0.055082 


0.055082 


0.0000 


0.925000 


0.925000 


0.0000 


-0.8 








0.0000 


0.900000 


0.900000 


0.0000 


-0.5 


1.743750 


1.743750 


0.0000 


0.750000 


0.750000 


0.0000 


-0.2 


2.700000 


2.700000 


0.0000 


0.600000 


0.600000 


0.0000 


0.1 


2.936250 


2.936250 


0.0000 


0.450000 


0.450000 


0.0000 


0.4 


2.520000 


2.520000 


0.0000 


0.300000 


0.300000 


0.0000 


0.7 


1.518750 


1.518750 


0.0000 


0.150000 


0.150000 


0.0000 


1 








0.0000 








0.0000 



Table 3b Comparison of velocity and temperature with Br = 0.05 , GR T = 5 , GR C = 5 , p = -5 and y* = -0.8 . 





Velocity 


Temperature 


y 


DTM 


PM 


Error 


DTM 


PM 


Error 


-l 








0.0000 


1.000000 


1.000000 


0.0000 


-0.95 


0.056795 


0.056380 


0.0004 


1.019848 


1.006548 


0.0133 


-0.9 


0.075886 


0.075222 


0.0007 


1.039612 


1.013017 


0.0266 


-0.85 


0.057042 


0.056460 


0.0006 


1.059364 


1.019475 


0.0399 


-0.8 








0.0000 


1.079032 


1.025854 


0.0532 


-0.5 


2.076642 


1.976262 


0.1004 


1.070430 


0.974183 


0.0962 


-0.2 


3.226545 


3.067590 


0.1590 


0.925078 


0.826975 


0.0981 


0.1 


3.511194 


3.337492 


0.1737 


0.748205 


0.658004 


0.0902 


0.4 


3.009255 


2.861363 


0.1479 


0.565508 


0.485024 


0.0805 


0.7 


1.804303 


1.717982 


0.0863 


0.341631 


0.283674 


0.0580 


1 








0.0000 








0.0000 


»c Comparison of velocity and temperature with Br = 0.09 , GR T = 5 , GR C = 5 , p = -5 and y* 




Velocity 


Temperature 


y 


DTM 


PM 


Error 


DTM 


PM 


Error 


-1 








0.0000 


1.000000 


1.000000 


0.0000 


-0.95 


0.061213 


0.057167 


0.0040 


1.160698 


1.031787 


0.1289 


-0.9 


0.082935 


0.076482 


0.0065 


1.320672 


1.063431 


0.2572 


-0.85 


0.063199 


0.057563 


0.0056 


1.480419 


1.095055 


0.3854 


-0.8 








0.0000 


1.639913 


1.126537 


0.5134 


-0.5 


3.134164 


2.162271 


0.9719 


2.084440 


1.153530 


0.9309 


-0.2 


4.901119 


3.361663 


1.5395 


1.958644 


1.008556 


0.9501 


0.1 


5.341092 


3.658485 


1.6826 


1.698429 


0.824406 


0.8740 


0.4 


4.567218 


3.134453 


1.4328 


1.413275 


0.633044 


0.7802 


0.7 


2.713644 


1.877367 


0.8363 


0.952144 


0.390614 


0.5615 


1 








0.0000 








0.0000 
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Table 4a Comparison of velocity and temperature with Br = o , GR T = 5 , GR C = 5 , p = -5 and >•* = 0.8 . 





Velocity 


Temperature 


y 


DTM 


PM 


Error 


DTM 


PM 


Error 


-l 








0.0000 


1.000000 


0.850000 


0.0000 


-0.7 


2.720194 


2.720194 


0.0000 


0.700000 


0.550000 


0.0000 


-0.4 


4.232842 


4.232842 


0.0000 


0.400000 


0.250000 


0.0000 


-0.1 


4.649777 


4.649777 


0.0000 


0.100000 


0.100000 


0.0000 


0.2 


4.052842 


4.052842 


0.0000 


0.075000 


0.050000 


0.0000 


0.5 


2.495194 


2.495194 


0.0000 


0.025000 





0.0000 


0.8 








0.0000 


1.000000 


0.850000 


0.0000 


0.85 


0.019844 


0.019844 


0.0000 


0.400000 


0.250000 


0.0000 


0.9 


0.026250 


0.026250 


0.0000 


0.100000 


0.100000 


0.0000 


0.95 


0.019531 


0.019531 


0.0000 


0.075000 


0.050000 


0.0000 


1 








0.0000 


0.025000 





0.0000 



Table 4b Comparison of velocity and temperature with Br = o.oi , GR T = 5 , GR C = 5 , p = -5 and > * = 0.8 . 





Velocity 


Temperature 


y 


DTM 


PM 


Error 


DTM 


PM 


Error 


-l 








0.0000 


1.000000 


1.000000 


0.0000 


-0.7 


2.825637 


2.816608 


0.0090 


0.924245 


0.918020 


0.0062 


-0.4 


4.412157 


4.396750 


0.0154 


0.797661 


0.789269 


0.0084 


-0.1 


4.859428 


4.841378 


0.0181 


0.658724 


0.649363 


0.0094 


0.2 


4.243934 


4.227453 


0.0165 


0.517827 


0.507653 


0.0102 


0.5 


2.615207 


2.604832 


0.0104 


0.362559 


0.352732 


0.0098 


0.8 








0.0000 


0.160579 


0.155260 


0.0053 


0.8 








0.0000 


0.160579 


0.155260 


0.0053 


0.85 


0.020506 


0.020448 


0.0001 


0.120437 


0.116447 


0.0040 


0.9 


0.027007 


0.026941 


0.0001 


0.080292 


0.077632 


0.0027 


0.95 


0.020005 


0.019963 


0.0000 


0.040147 


0.038817 


0.0013 


1 








0.0000 








0.0000 


4c Comparison of velocity and temperature with Br = 0.05 , GR T = 5 , GR C = 5 , p = -5 and y* 




Velocity 


Temperature 


y 


DTM 


PM 


Error 


DTM 


PM 


Error 


-1 








0.0000 


1.000000 


1.000000 


0.0000 


-0.7 


3.707230 


3.202265 


0.5050 


1.536854 


1.190102 


0.3468 


-0.4 


5.914596 


5.052384 


0.8622 


1.616169 


1.146343 


0.4698 


-0.1 


6.618204 


5.607785 


1.0104 


1.571084 


1.046814 


0.5243 


0.2 


5.848729 


4.925897 


0.9228 


1.507863 


0.938267 


0.5696 


0.5 


3.624481 


3.043386 


0.5811 


1.314728 


0.763661 


0.5511 


0.8 








0.0000 


0.674938 


0.376298 


0.2986 


0.8 








0.0000 


0.674938 


0.376298 


0.2986 


0.85 


0.026132 


0.022866 


0.0033 


0.506220 


0.282234 


0.2240 


0.9 


0.033437 


0.029704 


0.0037 


0.337487 


0.188161 


0.1493 


0.95 


0.024023 


0.021690 


0.0023 


0.168751 


0.094085 


0.0747 


1 








0.0000 








0.0000 
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Appendix 

1 1 

c, = --, c 2 = - ,c, 

2 2 



Sink (a y * ) + n Sink (a ) 



1 1 

-, C 4 = -, B l = 

2 2 Sin^ (a y *) Cosh (a ) + Sin/j (a J Cosh (a y *) 



n C o sh [a } - C o sh (a y * ) 



S in A (a y *)Cos/i(a ) + S in « (a ) C o s /? ( a y *) 



(p-Gi? r C 2 ) Gi^C, 
r. = , r, = 



Gfl fi, 



G* r C 3 



U (y * 2 -l) + r 2 (y * 3 +l) + r 3 (Cosh (a y *) - Cosh (a ))+ r 4 (Sinn (a y *) + Sinn (a ))) 



1 + y 



A 2 = Aj - r y + r 2 - r 3 Cosh (a ) + r 4 Sinii (a ) , A 3 



(i- y* 2 )+ r 6 (i- y* 3 ) 

(y*-l) 



A 4 = ~\ - r 5- r 6 



(2 A, 2 + r 4 2 a 2 - r*a 2 ) 2 A r ( 4 ^ + 6 A i r i) 3 r r. 



P 6 = " " . P 1 



3 12 
(2A,r 4 « ! -8r | r 3 a + 36r 2 r 4 ) 



Pa 



10 



( 2 A, r 3 a 2 - 8 r 4 a + 36 r 2 r 3 ) (4 ^ r 4 a - 24 r 2 r 3 ) (4 ^ r 3 a - 24 n, r 4 ) 

P 1 = ~ 5 ' PlO = ~ 2 ' Pi 1 = ~~ 2 



6 r 2 r 4 
Pl2 = ' Pl3 



6 r 2 r, A~ 

' 9 > = ~~T 

a 2 



2 A 3 r 5 



(2r/ + 3A 3 r 6 ) 



9 4 



3 r r 

J 5 6 



3r 



3 6 5 

I Pi — P 2 + P 5- P \ P + 5 P C g sh (2. cc ^) — p Siqh (2 ar ) -+- p C ogh (a ) \ 

^5 ' T = — 

I — p 9 Sinh (a ) ~ p l £osh(a)+p fifnh^aj+p C^o^sh (a ) — /? Siiph (a ) 

r 2 = ~ (?1 + <7 2 + 9 3 + 94 + <? 5 ) . 

T 3 = q i y * 2 + <? 2 >■ * 1 + 9 3 y * 4 + <? 4 y *' + 9 s >■ * 6 - p , y * 2 - p 2 >■ * 3 - p 3 y * 4 - p 4 y * 5 - p s y * 6 

— /? 6 C o i'/z ( 2 a y * ) — p Si'n/i (2« v *) - /? s C o sh ( a y * ) — p 9 S in h ( a y * ) 
~ P 10 y * C o sh (<z y 111 ) — p^y^Sin/jfary*) — Coi/i (a *) - p 13 y* 2 S/n/z(a;y*) 

T 4 = 2 # , y * + 3q 2 y * 2 + 4 3 y * 3 + 5 <jr 4 y * 4 + 6 5 y * 5 - 2 p t y * - 3 p 2 y *" - 4 /? 3 y * 3 -5 p t y * 4 - 6 /;> 5 y * 5 
- 2a 5in/i (2 a y *) - 2a Co.sft(2a_y*) - /? 8 a 5m/i(aj'*) — a p v C o sh (a y *) 
- p l0 ( y * a S in h (a y *) + Co sh [a >'*)) — p ,,(>'* a Co.s-ft(a>'*) + S/n/i(av*)) 
— ! , (2 y * C o .v /( ( a y * ) I « y * 7 SmA (ay*)) - p 13 (2y*Si'nA(ay*) + ay* 2 Coifi (a y *)) 

( y * r 4 + r, - r 2 - r 3 - r 4 ) (r, + r 2 + r 3 + r 4 (1 - y *)) 

G, = ,G, = , 



(-r, + r 2 + r 3 -r 4 (i+ y *)) 



^ GR T p, ^ = _G^ I _p 2 



g . = r - g, 



Gi? r G 2 



1 2 

G R r P 6 
4a 2 



20 



6 

GR T P 5 



3 



4a 2 



42 56 
(p 8 a 2 -2 P|1 a +6 Pl2 )GR I 
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(p 9 a 2 - 2 P l0 a + 6p 13 )Gi? r 



p l0 a -Ap n )GR 7 



p n a -4p 12 )GR T 



GR T P l2 

2 

a 

GR T q. 



GR T p n 



GR T G 3 



GR T q 3 



GR T <?. 

12 
GR T q 5 



20 30 42 56 

( /?, - R 2 + R } - R 4 + R f - R 6 + « 7 + R s Cosh (2 a) - R 9 Sinh (2 a ) + R l0 Cosh (a )"| 

R u Sinh (a) - R n Cosh (a ) + /? 13 Sinh (a) + R u C osh (a ) - R l5 Sinh (a ) 

f R t y * 2 +R 2 y * 3 +R 3 y ** + K 4 y * 5 + R 5 y *" + « 6 y * ? + « 7 y ** + /? 8 Co.f/i(2ay*) + R 9 5('n/i(2«y*) A | 



r 7 =1 +R io Cosh (a y*) + R n Sinh (a y *) + « 12 y * Co.s/z (« y *) + K,,y * Sinfc (ay*) 
+ i? 14 y * 2 Cos/i (a y *) + « 15 y * 2 Sinh (a y *) 



i + y 



1-y 



T 5 + G 5 , G s 



Issn 2250-3005 



|October| I2013 



Page 44 



